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^ I The ion phase-space dynamics in the Coulomb explosion of very large (~ 10® — 10^ atoms) deu- 

CIh' terium clusters can be tailored using two consecutive laser pulses with different intensities and an 

O ■ appropriate time delay. For suitable sets of laser parameters (intensities and delay), large-scale 

j>^, shock shells form during the explosion, thus highly increasing the probability of fusion reactions 

CL|' within the single exploding clusters. In order to analyze the ion dynamics and evaluate the in- 

tracluster reaction rate, a one-dimensional theory is used, which approximately accounts for the 
vQ , electron expulsion from the clusters. It is found that, for very large clusters (initial radius ~ 100 

nm), and optimal laser parameters, the intracluster fusion yield becomes comparable to the inter- 

(^ • cluster fusion yield. The validity of the results is confirmed with three-dimensional particle- in-cell 

^. 

^D . simulations. 
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I. INTRODUCTION 



Nowadays, CPA (chirped pulse amplification) laser technology jjj allows the table-top 
production of ultra-short (~ 10 — 100 fs) laser pulses, with peak power up to the Petawatt 
level J2| . When focusing such lasers to spot sizes of a few tens of micrometers, peak intensities 
up to 10^^ W/cm^ can be achieved, opening new research realms in the field of light-matter 
interaction. Among these, the interaction of ultra-intense lasers with clustered gases has 
become . cental .e.ea.cK .op.c BQQaflaflQQ. ..n,, o„.. to .he a.az..ly 
efficient coupling of ultra-intense laser radiation to clustered media; in experiments, nearly 
100% of the total laser energy has been observed to be deposited within a few millimeters 
propagation length J9|, |l^. This effective energy absorption results in various experimental 
evidences, such as bright x-ray emission [13|, [ij, |l5|, Il6 |. prod uction of highly-ionized matter 
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3, 



ISlllQll, generation of energetic electrons and ions [20, 
„ffusio„„eu.o„sflBfla. 
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and copious production 



Depending on the gas, the size of the clusters, and the laser features, a huge variety 
of physical scenarios is possible, from the slow, hydrodynamic expansion of quasi-neutral 
nanoplasmas to the violent Coulomb explosion of highly-charged ion clouds. In this paper, 
the attention is focused on the nonlinear ion dynamics in the Coulomb explosion of large 
(10^ — 10^ atoms) deuterium clusters irradiated with ultra-intense (10^*^ — 10^^ W/cm^) 
lasers. In the conditions considered here, all the atoms in a cluster are immediately stripped 
of their only electron, via field ionization, at the leading edge of the laser pulse, and the 
nanoplasma approximation [l7| can be adopted. The newly- formed free electrons can then 
be expelled, partially or totally, from the host cluster (outer ionization). If this process 
is extremely fast with respect to the explosion time scale, its dynamics may be neglected, 
but such an extreme scenario is met only with small clusters (less than 10000 atorns). For 
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26| and 



large clusters (10^ — 10^ atoms), the outer ionization dynamics is crucial 

it highly affects the explosion features. In the present work, we explore the possibility of 

using two sequential laser pulses ( "double pump" henceforth) to control the outer ionization 



process and to drive large-scale shock shells [2^|27| during the Coulomb explosion of large 
deuterium clusters. As a consequence of the shock shell formation, relative velocities appear 
between ions belonging to a single exploding cluster and intracluster DD fusion reactions 
may take place in the early phase of the explosion, long before the neighboring clusters start 



interacting with one another. It is found that, for optimal double-pump parameters, the 
intracluster fusion yield can become comparable with the intercluster fusion yield. Hence, 
for appropriate experimental conditions, a time-resolved burst of fusion neutrons should 
be detected before the usual bulk of fusion neutrons, thus providing a clear experimental 
evidence for the formation of shock shells on a nanometer scale J27|. 

In the present work, the laser-induced Coulomb explosion is modeled using an approx- 
imate ID theory, in which a prescribed electron dynamics, related to the pulse envelope 



function, is assumed [2^. The validity of the results is then confirmed with highly-realistic, 
three-dimensional (3D) particle-in-cell (PIC) 28j. The simulation framework adopted here 
is OSIRIS 2.0 29|, a state-of-the-art, massively-parallel, electromagnetic, fully-relativistic, 
3D PIC code. 

This paper is organized as follows. In Section |nj the formation of large-scale shock shells 
in Coulomb explosions is introduced. In Section ITTTl a ID model for the Coulomb explosion is 
outlined, with an emphasis on the inclusion of the effects caused by the electron dynamics. 
In Section IIV| the ID model is employed to study the intracluster fusion reactions. The 
fusion yield from intracluster and intercluster reactions are compared in Section and 
three-dimensional PIC simulations of the double-pump technique are presented in Section 
fVll Finally, the conclusions are stated. 



II. LARGE-SCALE SHOCK SHELLS IN COULOMB EXPLOSIONS 

The possible formation of shocks J30| during the Coulomb explosion of an ideal pure-ion 

I 
sphere has been recently demonstrated theoretically \T^. In particular, it has been shown 

that every radial nonuniformity in the initial ion density profile leads to the formation of 
a multi-branch phase-space structure, which is named a shock shell, accompanied by the 
appearance of one or more peaks in the ion density profile. Such structures form when the 
initial ion density is nonuniform and a maximum in the radial velocity profile soon forms at 
some radius within the distribution, leading to a characteristic overtaking process between 
the ions. 

One interesting consequence of the shock formation is the possible occurrence of ener- 
getic ion-ion collisions within a single expanding cluster. For instance, the appearance of 
pronounced shock shells in the explosion of very large (10^ — 10^ atoms) deuterium clusters 



could cause an enhancement in the probabihty of intracluster fusion reactions. According 
to the theory, in order to obtain a large-scale shock shell, the nonuniformity in the initial 
density profile must be large-scale as well. In fact, if the initial density profile differs just 
slightly from a uniform, step-like profile, the corresponding shock shell is doomed to be 
small-scale and cannot involve significant relative velocities and collision energies. Results 
from 3D numerical simulations of the interaction of an ultra-intense laser with large deu- 
terium clusters in realistic conditions showed that shocks may form naturally even when 
starting from uniform, step-like density profiles J2J] because the electron dynamics in the 
laser field smooths out the ion density profile in the early stages of the interaction. However, 
the simulations also showed these spontaneously-occurring shocks are small-scale, so that no 
appreciable relative velocities appear within the cluster, which means that no shock-driven, 
intracluster fusion events occur. 

At first sight, inducing such reactions via generation of large-scale shock shells in ex- 
ploding deuterium clusters would then appear unfeasible, as it would require the ability to 
tailor the initial radial density profile of roughly spherical objects on the nanometer scale. 
Surprisingly, on the contrary, such an apparently-difficult goal can be achieved quite simply 
by using a suitable sequence of laser pulses with different intensities [2.4]. In fact, if the first 
pulse is relatively weak, only a small number of electrons are stripped off the host cluster, and 
a slow expansion takes place, driven by both Coulomb repulsion and hydrodynamic pressure 
of the hot electrons jl2|. As the expansion goes on, the ion density in the outer region of 
the cluster decreases, while the cluster core remains dense and approximately neutral: in 
this way, a smoothly-decreasing plasma density profile is naturally formed. Afterwards, if 
the electrons are suddenly swept away from the cluster core by a second, extremely-intense 
pulse, the inner ions immediately feel a much higher repulsive force than the outer ions do. 
The cluster core explodes abruptly and the fastest inner ions overrun the slowly-expanding 
outer ions, leading to the formation of a large-scale shock shell. Such a situation is well 
depicted in Fig. ^ where colored spheres indicate position and energy of a random sample 
of ions, after the passage of the second pulse (the figure has been built from the results of 
the 3D PIC simulation discussed in Section IVH) . 

The double-pump technique not only provides an effective method to induce large-scale 
shock shells in the Coulomb explosion of very large clusters, but it also provides the ability 
to control the phase-space structure of the exploding clusters, by varying the delay between 



the two pulses, At, and the peak intensity of the first pulse, /i. The dependence of the 
explosion features on these parameters is analyzed in detail in the following sections, where 
we investigate the possibility of inducing intracluster, shock driven fusion reactions in double- 
pump experiments. 

Alternatively, large-scale shock shells can also be driven using clusters constituted by 
more than one species, as long as the ion charge-to-mass density ratio decreases radially. 
Such possibility will be explored in a future publication J31|. 

III. ID MODEL FOR COULOMB EXPLOSION 

A proper tuning of the double-pump parameters requires a deep knowledge of their influ- 
ence on the explosion dynamics, the shock shell formation and evolution, and, consequently, 
the intracluster fusion yield. Since a direct parameter scan via either 2D or 3D PIC sim- 
ulations would have been impracticable, because it is too computationally demanding, we 
have developed a simple ID theoretical model that allowed us to investigate the effects of 
the key double-pump parameters, namely the time delay between the pulses. At, and the 
peak intensity of the first pulse, Ji. 

Our model, which is an extended version of the standard model for the Coulomb explo- 
sion of a spherical pure-ion cluster (cf. J27|), is suitable for both spherical and cylindrical 
symmetry and takes into account the effect of the electron population in the expanding 
cluster. The electron population follows a prescribed dynamics, determined uniquely by the 
laser pulse features. In this framework, the acceleration of an ion at a given time r and 
radial position r can be written, in dimensionless units, as 

£r _ g(r)-gei(r,r) .. 

where 7 accounts for the geometry (7 = 2 for spherical geometry, 7 = 1 for cylindrical 
geometry), and where mass is normalized to m (ion mass), length to Rq (initial radius of the 



cluster), charge to the elementary charge e, and time to to = J [mWo^ ) / (e^A^o); ^0 also 
represents the time scale for the explosion of a pure-ion spherical/cylindrical distribution, 
being A"o the total number of ions (per unit length in the cylindrical case). The quantity 
Q{r) is the ion charge within a sphere/cylinder of radius r, while (5ei(^, t) describes the 
absolute value of the electron charge still present within the same sphere/cylinder at time r. 



Using the Cluster Barrier Suppression Ionization (CBSI) model described in J32], Qc\{r,T) 

can be determined from the instantaneous value of the laser-field envelope function, Ei{t). 

According to the CBSI theory, an electron is expelled from the cluster boundary to infinity 

whenever 

Ei{t) l-Q4R,r) 

where R is the cluster radius at time r (the \/2 factor accounts for the periodicity of the 
laser electric field). To determine Qe\{fiT) we need a further assumption on the radial 
profile of the electron distribution. The simplest choice is to assume that, at each time r, 
the remaining electrons fully neutralize the cluster core. The permanence of a core of cold 



recently explained theoretically in 

Qei(r,r) : 
being 



M 



electrons within the cluster is clearly visible in PIC simulations 

34| . Thus, Qe\{f-iT) takes the form 
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the radius of the electron sphere/cylinder at time r. Under these assumptions, Eq. (^ can 
be written as 

■d^2 = \ Q{r)-[l-R'<Ei{T)/V2\ 



(5) 
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r > i?el(^ 



providing a simple model for studying the Coulomb explosion of large clusters driven by a 
general sequence of laser pulses, having different intensities and envelopes. This model can 
be used to analyze scenarios involving the formation of large-scale shock shells, allowing the 
control of the shock features through tuning of the pulse parameters. 

As long as the initial ion density profile po{ro) is known, Eq. (jH)) can be readily integrated 
numerically by following the trajectory r(ro, r) of a finite set of ions with different initial 
position tq. Details of the numerical model will be described elsewhere. 

IV. INTRACLUSTER FUSION REACTIONS 



Large-scale shock shells are characterized by a well-defined multi 
V — r phase space, most frequently a three-branch structure J2J, |27| 



)ranch structure in the 
as the one in Fig. 12 



which refers to the Coulomb explosion of a pure-ion sphere with a nonuniform radial density 



profile [27[. As the explosion goes on, the upper branch overlaps the lower branches: the 
shock shell widens radially, narrowing its velocity spread, and, meanwhile, the ion density 
on each branch decreases. Therefore, one can reasonably expect the probability of nuclear 
reactions between ions belonging to different branches to be higher in the early stages after 
the shock shell formation and to decrease rapidly at advanced times. 

At each radius r and time r, the number of reactions per unit time and unit volume, 7?., 
is given (in dimensionless units) by 

n=l f f f {r,v) f {r,v') a {\v - v'\) \v ~ v'\dvdv' , (6) 

2 Jv Jv' 

where / (r, v) is the ID distribution function for the ions and a the normalized DD fusion 
cross section {a = ct/Rq, with a in cm^). Outside the shock shell, this integral vanishes 
because there are no relative velocities; within the shock shell, it simply reduces to the sum 
over the three branches of the phase space profile (identified in color in Fig. 12)) : 

1 ^ 

n= -^Piir)pj{r)a{\vi-Vj\)\vi-Vj\ (7) 

where Pi(r) = 1/ {2'~fTTr"') dQi/dr is the ion density on the ith branch. The intracluster 
reaction rate, R, is then 

R = Nl I ^"^ n 277rr^dr, (8) 

where r^h and i?sh represent the shock shell boundaries. The number of reactions per cluster, 
A/", is given by 

/•oo 

M= Rdr (9) 

where Tsh is the shock formation time. Once the phase space history of each ion of the initial 
distribution is known, from the solution of Eq. (jj)). A/" can be evaluated numerically through 
Eqs. ©-0. 

The ID theory here described provides a useful framework to perform parametric studies 
and investigate the influence of Ji and At on the total number of reactions per cluster. A/", 
seeking the combination of parameters that maximizes it: a maximum in A/" is expected to 
show up when the first pulse is intense enough to drive the first expansion but not so intense 
to expel too many electrons from the cluster. The delay of the second pulse needs to be 



long enough to allow the formation of a decreasing density profile but not so long to let the 
outer ions expand to large r, far from the cluster core. 

When analyzing the results of the parametric studies, the approximations introduced in 
the ID theory must be kept in mind. First, as hydrodynamic effects are neglected, the 
effect of the first pulse is underestimated when dealing with low values of Ji, meaning that, 
in reality, the first expansion is faster than what is predicted by Eq. (0) and the optimal 
value of Ji is actually lower than expected. Furthermore, the ID theory cannot include 
laser polarization effects which may affect the dynamics by causing an asymmetry in the 
explosion. Finally, the ID theory with CBSI model is expected to break down with enormous 
clusters {Rq > 500 nm - 1 /im), having initial radius similar to the laser central wavelength, 
meaning that the predicted optimal double-pump parameters could be unreliable in such 
cases. 

We have tested the validity of the reduced model by comparison with a series of results 
from single-pump 2D PIC simulations in various configuration of cluster size and laser peak 
intensity. Figures Ek andEb show comparisons between lineouts (along the laser propagation 
(x) and polarization (y) directions) of the ion density distribution, taken from the simulation 
results, and the density profile predicted by the ID theory. Both plots refer to the interaction 
of a circular, rodlike cluster (initial radius Rq = 32 nm, particle density no = 4.56 x 
10^^ cm~^) with a laser pulse having central wavelength Aq = 820 nm, and approximately 
Gaussian envelope with rise time trise = 35 fs. The peak intensities are /a = 4 x 10^^ W/cm^ 
(Fig. EK) and Jb = 1.6 x 10^^ W/cm^ (Fig. Eb)- As expected, the results from the ID 
model are quantitatively more accurate in the high-intensity case, where Coulomb forces 
are dominant, than in the low- intensity case, where hydrodynamic effects and polarization 
effects are relevant, and the cluster expands more rapidly than what is predicted by the ID 
theory. For these reasons, when exploring the possibility of intracluster fusion reactions, we 
resorted to the ID model to seek a good combination of double-pump parameters, which we 
then adjusted and employed to perform accurate 3D PIC simulations presented in Section 

ED 

The model also provides information on how the initial cluster size Rq affects A/", allowing 
one to perform PIC simulations using clusters with initial radius Rq ~ 10 nm and then to 
extrapolate the results to the case of extremely large clusters having initial radius Rq ~ 100 
nm, without performing new, and extremely large, simulations. 

8 



Here, we present results from parametric studies, with respect to Ji and At, for several 
cluster sizes, with initial radii, Rq, in the range 16 — 200 nm. For each Rq, we consider a 
spherical cluster of atomic deuterium (with uniform step-like density profile, no = 4.56 x 10^^ 
cm~^) hit by a pulse sequence in which a weak laser pulse (peak intensity Ii (variable), central 
wavelength Ao,i = 820 nm, and approximately Gaussian envelope with rise time trise,i ~ 35 fs) 
is followed by an ultra-intense pulse (peak intensity I2 ^ h, central wavelength Ao,2 = 820 
nm, and approximately-Gaussian envelope with rise time trise,2 ~ 20 fs) with time delay 
At variable in the range 70 — 500 fs. In all cases, the peak intensity of the second pulse 
is assumed to be high enough to expel all the electrons from the cluster core and drive a 
sudden Coulomb explosion. Figure HI shows the total number of reactions per cluster. A/", as 
a function of /i and At, for two representative cases: Rq = 32 nm (Fig. 0^) and Rq = 100 
nm (Fig. HJo). With Rq = 32 nm. A/" assumes its maximum value, A/^nax = 4.85 x 10~^ 
reactions, for h = 8.6 x lO^^^ W/cm^, At = 236 fs. With Rq = 100 nm, A^max = 5.58 x lO^^ 
reactions for Ii = 1.4 x 10^^ W/cm^, At = 139 fs. Hence, approximately trebling the cluster 
size (from Rq = 32 nm to Rq = 100 nm) results in gaining three order of magnitudes in A/", 
with a first pulse sixteen times as intense and a much shorter delay. These drastic changes 
are partly due to the variation of the DD fusion cross section, a, with the collision energy. In 
fact, as follows from Eqs. ©-(0), if cr were constant, A/max would be proportional to Rq, the 
optimal intensity for the first pulse would scale as R^, while the optimal delay would stay 
the same. In that case, the two plots in Fig. HI would have the same shape and A^ax would 
increase by a factor less than 100. The dependence of the optimal combination of Ji and 
At on the initial cluster size is depicted in Fig. EJ along with the corresponding variation of 
ACax- As Rq increases, the optimal intensity of the first pulse increases as well, while the 
optimal delay decreases and appears to saturate towards the value At = 125 fs. 

For both cases showed in Fig. |3J we report the evolution of the v — r phase space profile 
(starting right after the second pulse reaches the expanding cluster causing the formation 
of the shock shell), along with the time history of the reaction rate R, in Figs. IHlandlZl re- 
spectively. As one can see, R exhibits a sharp peak immediately after the shock formation, 
which occurs at time ti, causing the probability of intracluster fusion reactions to retain 
appreciable values only for a time interval of a few tens of fs. Intracluster fusion reactions 
represent a much faster, and profoundly different, phenomenon than intercluster reactions 



(whose typical time scale is ~ 10 — 100 ps 



35|): they occur abruptly in the single clusters 



when the particular phase-space dynamics origins a tiny, high-density, and short-hved reac- 
tion volume (the shock shell). On the contrary, intercluster reactions occur later in the big, 
long-lived reaction volume which is the whole plasma filament created by the laser pulses. 
An estimate of the respective contributions of the two, distinct processes to the total fusion 
yield is presented below. 

V. INTRACLUSTER AND INTERCLUSTER FUSION YIELDS 

In the last years, various experiments revealed the occurrence of nuclear fusion reactions 
in clustered media irradiated by ultra- intense fs lasers y, H ly, |7| . Most theoretical models 



developed to explain the experimental results showed that fusion reactions arise 
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33, 



primarily 



35[, though 



from collisions between fast ions belonging to different clusters 
the contribution from collisions of fast ions with colder ions and atoms outside the plasma 
man,ent c.ea.ed .. t.e la., .a. aI.o been .ve.,ated H. He.e. we a„a„.e .Ke .ole of 
intracluster, shock-driven fusion reactions in double-pump scenarios, for different cluster 
sizes, and we compare the intracluster fusion yield with the intercluster fusion yield. 

For the calculation of the intercluster fusion yield, Yic, we refer to the simple model 
presented in J35|, where Yic is evaluated as 

>^ic = ^n^ Wv) VrTa (10) 

being n the average atomic density inside the reaction volume K, Tj the plasma disassembly 
time (typical time for the expansion of the plasma in the reaction volume), and 

D /»oo/'00/'7r 

<^^) = W P (S,) P [S,) a{S^^n) X 

■^1^0 Jo Jo Jo 

X {2£con/mY''^ sm{a) d£id£2da (11) 

where £^coii = Si + £2 — '^ (^1^2) cos(a) is the binary collision energy for particles with 
kinetic energy Si = vl/2, S2 = vl/2, and collision angle a. The quantity P (S) is the 
(dimensionless) energy distribution of a single exploding cluster, calculated as 

P{S) = 11 f (r', v') S (v''^l2 - S\ 47rr'^dr'dw' = 
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1 dQ 

V\ dv 



' ^ 1)- rli) 



(12) 
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being {v^} the N^ zeros of v'^/2 — S. The reaction volume is assumed to be a cyhnder with 
radius Rr = 100 fun and height Hj. = 2 mm, and the plasma disassembly time is estimated 



as Td ~ 7r/2r/ (2t;max) |35|. 



In the previous section, we analyzed the explosion dynamics of a single cluster, neglecting 
its interaction with the neighboring clusters, and calculated the probability of intracluster 
reactions. Now, and in order to compare intracluster and intercluster fusion yields, we must 
consider the whole spatial distribution of clusters. For the sake of simplicity, we consider the 
ideal situation of a reaction volume containing a fixed number of deuterium atoms grouped 
in A^ci spherical clusters of equal size -Rq, with a mean separation distance d, distributed on 
a regular square lattice. Given the average atomic density n, the total number of clusters is 
A'^ci = nRlHru/No and the intercluster distance is d = Rq [47rno/ (3n)] ' . The intracluster 
fusion yield. Fie, can be evaluated as 

Fie = iVeiAT (13) 



provided that the ratio 6 = d/ Rq (which depends only on the packing fraction n/riQ 33] ) be 
sufficiently high (5 > 20 — 25) and, consequently, most intracluster reactions occur before 
the cluster cores start interacting with one another. In the opposite case (small values of 
5: (5 < 10 — 15), almost no intracluster reactions can occur before the cluster cores start 
interacting with one another. This is explained by noticing that, typically, the intracluster 
reaction rate, R, peaks when the radius of the outer boundary of the shock shell, i?sh) 
is i?sh ~ 4 — 5-Ro) and decreases below 1/10 the peak value when _Rsh ~ 10 — 15-Ro (see 
Figs, ini IZ|)- If ^^ = 10^^ cm~^, a typical experimental value, and with the cluster density 
riQ = 4.56 X 10^^ cm~^ here considered, one obtains d/Ro — 27. In such a low packing fraction 
case. Fie can be evaluated through Eq. (jl3|) and a comparison with Fic is readily carried out. 
Figure |H1 reports the value of Fjc and Yjc for different cluster sizes: for each value of i?o "we 
sought the optimal combination of double-pump parameters (see Fig. EI) and then calculated 
the corresponding fusion yields. The contribution of the intracluster reactions grows rapidly 
with Rq, meaning that the growth of the shock shells and the increase of the fusion cross 
section prevail against the decrease of the number of clusters in the reaction volume. On the 
contrary, the intercluster fusion yield, though keeping a high value, decreases at very high 
Ro because collision energies beyond the one that maximizes a start to appear. This is well 
illustrated in Fig. El where the maximum intracluster and intercluster collision energies are 

11 



compared for different cluster sizes, along with a (£^coii)- For initial cluster radii > 70 — 80 
nm, the maximum collision energy for intercluster reactions lies to the right of the peak in 
o (£^coii); while the maximum collision energy for intracluster reactions lies to the left of it 
even for initial cluster radii as high as 200 nm: this explains the results of Fig. |H1 

From our discussion, we conclude that, in principle, a double pump experiment with 
optimized pulse parameters, very-large clusters, and low packing fraction, should provide a 
clear signature for the occurrence of intracluster, shock-driven fusion reactions, in the form 
of a time-resolved burst of fusion neutrons which anticipates the bulk of fusion neutrons 
produced via both intercluster reactions within the plasma filament and ion-ion/ion-atom 
collisions outside the plasma filament. 

VI. THREE DIMENSIONAL PIC SIMULATIONS 

In order to check the validity of the analysis presented above, and to get deeper physical 
insights, we have performed 3D PIC simulations of the laser-cluster interaction in a double- 
pump case, treating self-consistently the dynamics of electrons and ions in the laser field, the 
outer ionization dynamics, and the full dynamics of both the slow expansion induced by the 
first pulse and the sudden explosion driven by the second pulse. In PIC simulations, a set 
of computational particles is moved under the action of their self-consistent electromagnetic 
field and any externally applied field: this is done by first depositing the current density on 
a spatial grid, then solving Maxwell's equations on the same grid and computing the force 
accelerating each particle, by interpolation of the field values on the position of the point 
particle. 

We first consider the simulation of the irradiation of a cluster having radius i?o = 32 nm 
and density no = 4.56 x 10^^ cm~^ with a pulse sequence whose parameters are the optimal 
ones calculated in Section H^ (Ji = 8.6 x 10^^ W/cm^, At = 236 fs), except for the peak 
intensity of the first laser, which has been lowered to Ji :i^ 2 x 10^^ W/cm^, to compensate 
for the underestimation of the expansion velocity in the ID model. The second pulse, having 
peak intensity I^ — 1-3 x 10^^ W/cm^ and a shorter pulse duration (trise = 20 fs), hits the 
cluster with time delay At = 236 fs (same as the optimal value calculated above). Single- 
pulse 3D simulations had already shown that peak intensities lower than I2 suffice to expel 
all cluster electrons before the peak of the pulse reaches the cluster J2J|, even though the 
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electron dynamics cannot be assumed as instantaneous on the explosion time scale. Both 
pulses are linearly polarized along the z direction, propagate in the i, direction, and their 
envelope is approximately Gaussian. The simulation box is cubic, with side Lbox = 1 /^m, 
discretized in a 420 x 420 x 420 uniform spatial grid and the number of particles per species 
is 6.4 X 10^, a value close to the actual number of atoms, for the configuration described. 

We concentrate our analysis on the ion dynamics in the second part of the simulation, 
when the more intense laser interacts with the slowly-expanding cluster, driving its Coulomb 
explosion: right before the interaction with the second pulse begins, the ion density and 
phase-space profiles appear as in Fig. E3 where the density profile is decreasing from the 
center toward the periphery, but the expansion is clearly asymmetric, being much faster 
along the laser polarization direction, as also testified by the phase space lineouts. As a 
comparison, the density and velocity profiles obtained with the ID model are also plotted 
(in the optimal case I\ = 8.6 x 10^^ W/cm^, At = 236 fs). In the PIC simulation, the cluster 
has expanded slightly more, despite the lower value of Ji (meaning hydrodynamics effects 
are indeed very relevant for the configuration considered), and the density profile is different 
from the theoretical one, being the cluster core less dense. These differences, which are due 
the limitations of the ID model, where a crude approximation on the electron dynamics is 
adopted, do not prevent the formation of a large-scale shock shell. Yet, they appear to affect 
the Coulomb explosion dynamics, especially in limiting the maximum energy acquired by 
the inner ions (120 keV instead of the 210 keV predicted by the ID theory). This is clearly 
visible in Fig. ^J where the ion phase-space history is shown. Again, two curves along the 
y and z directions are plotted and compared with the theoretical curves. The core explosion 
is much more symmetric than the first slow expansion, since the electrons are quite rapidly 
expelled from the cluster and their dynamics has a smaller infiuence on the ion dynamics 
than during the first expansion. The explosion predicted by the ID theory is more violent 
for a variety of reasons: first, when the second pulse hits the plasma and the explosion 
starts, the Coulomb energy stored in the cluster core, where much more charge is packed 
(cf. Fig. ITn|l . is higher (which also explains why the dynamics of the outer ions, the ones 
forming the lower branch of the phase space profile, resembles the numerical results more 
than the dynamics of the inner ions). Furthermore, in the PIC simulations, neutralization 
by the hot electrons expelled from the cluster, but remaining in the computational domain 
(where the total net charge is zero), also play a role, as well as 3D effects, anisotropies, 
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boundary effects arising when the ion cloud gets as big as the simulation box, and possible 
propagation effects in the underdense expanding plasma surrounding the dense cluster core. 
The simulations results show that the optimal double-pump configuration predicted by 
the ID theory represents a good approximation for inducing the formation of large-shock 
shells, capable of driving intracluster nuclear reactions in realistic cases, even though all the 
effects mentioned above can highly affect the expansion/explosion dynamics. Their influence 
on the reaction yield needs investigation, for instance through numerical calculation of the 
reaction rates during the PIC simulations, to be presented in a future publication. However, 
a first estimate based on Eq. ((Zj) suggests that differences in the ion dynamics like those of 
Fig. ^3 where PIC calculations predict a shock shell with velocity amplitude reduced by 
a factor of two with respect to the ID model, should result in intracluster reaction yields 
reduced by a factor of the same order (the actual reduction depends on the variations of a 
for the conditions considered). In the case analyzed here, results also seem to suggest that a 
more pronounced shock shell, with dynamics more akin to the theoretical model, would be 
obtained with a weaker first pulse and/or a shorter time delay, and with a stronger second 
pulse. These trends should hold with larger clusters [Rq ~ 100 — 200 nm) too, provided that 
conditions lie in the range of validity of the ID theory. 

VII. CONCLUSIONS 

The explosion dynamics of large deuterium clusters irradiated by sequential laser pulses 
having different intensities has been investigated, focusing on the formation and evolution of 
large-scale shock shells in the ion phase space, which lead to the occurrence of intracluster 
fusion reactions. The effect of the double-pump parameters (delay of the second pulse and 
intensity of the first pulse) on the number of intracluster fusion reactions has been analyzed 
using a ID model which approximately accounts for the outer ionization dynamics. After the 
optimal double-pump configuration has been found, for different cluster sizes, the optimal 
intracluster fusion yield has been calculated and compared with the intercluster yield, finding 
that intracluster reactions become important with very large clusters (radius ~ 100 nm). 
The optimal double-pump parameters obtained with the ID model have then been used to 
perform three dimensional PIC simulations, whose results confirmed the formation of well- 
pronounced shock shells with high relative velocities inside the single exploding clusters. 
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FIG. 1: (Color online) Ion distribution in configuration space at time t = 350 fs for the double- 
pump case described in Section IVTl The colored spheres indicate position and energy of a random 
sample of ~ 1 x 10^ (out of ~ 6.4 x 10^) particles. Color is proportional to energy, the lightest 
spheres corresponding to the maximum ion energy, i?max = 120 keV. 
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FIG. 2: Typical three-branched shock shell in phase space (thick, gray lines mark the different 
branches) |2J, |27| for the Coulomb explosion of a pure-ion sphere with nonuniform radial density 
profile. Units are normalized as in Section ITTll 
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FIG. 3: Ion number density when the cluster radius is ~ 5Ro: (a) low- intensity case, and (b) 
high-intensity case. Thick gray lines represent lineouts in the x (dark) and y (light) directions. 
Thin black lines refer to the solution obtained from the ID theoretical model. 



19 




FIG. 4: (Color online) Total number of reactions per cluster, AA, as a function of the peak intensity 
of the first laser pulse, /i, and the time delay between the pulses, At, for (a) Rq = 32 nm and (b) 
Ro = 100 nm. 
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FIG. 5: (a) Maximum value of the total number of reactions per cluster, A/^naxj and (b) optimal 
values for the double pump parameters /i (gray squares) and At (black circles), for different values 
of the initial cluster radius, Rq. 
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FIG. 6: (a) Phase space profile at times ti = 255 fs, t2 = 305 fs, t^ = 355 fs, and ^4 = 405 fs; (b) 
time history of the reaction rate, R, for Rq = 32 nm, and with the optimal combination of double 
pump parameters: /i = 8.6 x 10"*^^ W/cm^, At = 236 fs. 
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FIG. 7: (a) Phase space profile at times ti = 155 fs, t2 = 205 fs, ts = 255 fs, and ^4 = 305 fs; (b) 
time history of the reaction rate, R, for Rq = 100 nm, and with the optimal combination of double 
pump parameters: /i = 1.4 x 10"*^^. W/cm^, At = 139 fs. 
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FIG. 8: Intercluster fusion yield, l^c (black circles), intracluster fusion yield, l^c (black squares), 
and percentage value of l^c/^c (gi'ay bullets) for different values of the initial cluster radius, Rq. 
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FIG. 9: Cross section for the deuterium-deuterium fusion reaction, a, as a function of the colh- 
sion energy (full, black line); gray squares and gray circles indicate, respectively, the maximum 
intercluster and intracluster collision energy for different values of the initial cluster radius, Rq 
(reported on the right, gray axis). 
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FIG. 10: (a) Ion number density, and (b) phase space profile at times t = At = 236 fs. In (a), grey 
thick lines represent lineouts in the y (dark) and z (light) directions. In (b), grey points mark the 
position in the v — r phase space for those particles contained in a solid angle AQ ~ 0.1 sr around 
the y (dark) and z (light) directions. Thin black lines always refer to the solution obtained from 
the ID theoretical model. 
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FIG. 11: Phase space profile at times ti = 270 fs, t2 = 315 fs, t^ = 350 fs. The gray points mark 
the position in the v — r phase space for those particles contained in a solid angle AO, ~ 0.1 sr 
around the y (dark) and z (light) directions. The black lines refer to the solution obtained from 
the ID theoretical model. 
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